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Abstract. We study two-dimensional quantum Gaussian packets with a fixed value 
of mean angular momentum. This value is the sum of two independent parts: the 
‘external’ momentum related to the motion of the packet center and the ‘internal’ 
momentum due to quantum fluctuations. The packets minimizing the mean energy 
of an isotropic oscillator with the fixed mean angular momentum are found. They 
exist for ‘co-rotating’ external and internal motions, and they have nonzero correlation 
coefficients between coordinates and momenta, together with some (moderate) amount 
of quadrature squeezing. Variances of angular momentum and energy are calculated, 
too. Differences in the behavior of ‘co-rotating’ and ‘anti-rotating’ packets are shown. 
The time evolution of rotating Gaussian packets is analyzed, including the cases of a 
charge in a homogeneous magnetic field and a free particle. In the latter case, the 
effect of initial shrinking of packets with big enough coordinate-momentum correlation 
coefficients (followed by the well known expansion) is discovered. This happens due to 
a competition of ‘focusing’ and ‘de-focusing’ in the orthogonal directions. 
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1. Introduction 


Quantum Gaussian wave packets were considered at the dawn of quantum mechanics 
by Schrodinger [I], Kennard [2| and Darwin |3|. Since that time, different properties of 
such packets were studied by many authors. We are interested here in two-dimensional 
packets. They were discussed in detail, e.g., in original papers, books and reviews (4^9] , 
where the main emphasis was made on the squeezing properties, or in more recent 


papers 10 -16 , where the problems of quantum information (in particular entanglement 


properties ) were the focus of studies. 

The aim of the present paper is to study properties of Gaussian packets possessing 
nonzero mean values of the angular momentum operator. Different special cases of such 
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packets were considered by many authors for a long time. In particular, packets whose 
centers move along some circles arise naturally in the description of quantum charged 
particles placed in a uniform magnetic field. An example of such packets was given for 
the first time by Darwin [3|. Similar ‘circulating wave packets’ were considered later 
by Husirni [l7|. These packets were generalized and interpreted as coherent states of a 


charged particle (or an oscillator) in a magnetic held by Malkin and Man’ko 18,19 and 
other authors |20,2l]. Using such packets one can simplify substantially the analysis 
of many physical effects, in particular the Landau diamagnetism 22 . The Gaussian 


packets in a magnetic held were studied also in 23,24 . Two-dimensional coherent states 


in rotated frames were shown to be useful for the description of rotational properties 
of nuclei in the framework of the so called ‘cranked oscillator’ model 25,26 . More 


recently, rotated Gaussian packets or Gaussian packets in rotating frames were studied 
Mathematically equivalent ‘twisted’ paraxial optical beams with nonzero 


m 


27-31 


orbital angular momentum were studied, e.g., in 32-35 


During the past two decades, the so called ‘vortex states’ of light beams with 

In parallel, many 


nonzero orbital angular momentum received much attention 36-41 


authors studied vortex states of quantum massive particles 42 -48 . Recently such states 


were created in different experiments 49 -52 . It was pointed out in the cited papers that 


matter wave beams with orbital angular momentum can hnd applications in different 
areas, such as condensed-matter spectroscopy, electron microscopy and particle physics. 
In particular, they can be used in the study of magnetic properties of materials and for 
manipulating nanoparticles. Therefore a study of properties of rotating packets that 
were not considered earlier seems to be an interesting and timely task. 

The new question answered in this paper is as follows: what is the minimal mean 
energy of quantum packets with a fixed value of mean angular momentuml This question 
has sense, e.g., for a free isotropic oscillator or an isotropic oscillator placed in a 
homogeneous magnetic field (with a free particle moving in the magnetic field as a 
special case). 

Solutions to the stationary Schrodinger equation for the two-dimensional harmonic 
oscillator of mass M and frequency c a in the polar coordinates are known since the very 
beginning of quantum mechanics [53] : 

-0 n r m(r , (p ) = iVr |m| <f> (~n r ; \m\ + 1; fir 2 ) exp (-fir 2 /2 + ivrup) , fi = Moj/h. (1) 

Here m = 0, ±1, ±2,..., n r = 0,1, 2,... and <h(a; c; z) is the confluent hypergeometric 
function. Function ([Tj) can be written also in terms of the generalized Laguerre 
polynomials. Using the definition |54|[55] 


<L(-n- a■ z) = + iG)( z ) 

[ n ' a ' Z) T(n + « +1) " 1 j 

we obtain the following expression for normalized solutions (frequently called, especially 
in the optics applications, as Laguerre-Gauss functions): 


fpn r m(r,(p) = 


fin r 


(fir 2 )' 1 " 1 ^ 2 (fir 2 ) exp (—fir 2 /2 + inup) . 


it ( n r + |m|)! 


( 2 ) 
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From the well known form of energy spectrum 

E nrm = hw (1 + \m\ + 2 n r ) (3) 

it is clear that the minimal energy eigenvalue for a fixed eigenvalue of the angular 
momentum hm equals E m i n (rn ) — Huj (1 + |m|). But what is the minimal mean value 
of energy for superpositions of energy eigenstates with a fixed mean value of angular 
momentum hC and an arbitrary value of £? Taking into account the orthogonality of 
functions ^n r m one can see that the answer is 

(E) min (£) = hcu( 1 + |£|). (4) 

This value is achieved, for example, in superpositions of states with zero value of radial 
quantum number n r and the same signs of quantum numbers m: 

'Ipmin ^ ^ t'rrpiuhn, • ^ ' | C m \ 1; ^ ' 7TI | C m | llC, ) ^ ' \l7l\ \ C rn h|T|. (o) 

mm m m 

Obviously, the number of possible superpositions of this kind is infinite. We wish to 
know, if there exist Gaussian packets satisfying equation Q? The answer to this 
question is positive, but it is not obvious or trivial, as we show in the subsequent 
sections. In particular, the result depends crucially on the mutual directions of the 
‘external’ rotation (related to the motion of the center of packet) and the ‘internal’ one 
(related to the evolution of quantum fluctuations and the direction of rotation of the 
ellipse of constant probability density). Note that the intensity (probability density) of 
the eigenfunctions |^ 0m | 2 with m^O equals zero at the center. In contrast, Gaussian 
packets with nonzero mean angular momentum have the maximal probability density 
at the center. Moreover, it will be shown in section [3] that the minimizing packets 
maintain their shape and size in the coordinate space during the time evolution (for a 
nonzero oscillator potential or nonzero magnetic field), rotating like a rigid body, but 
being squeezed at the same time. This is a generalization of Schrodinger’s packets [I] 
that gave rise eventually to the concept of coherent states. Therefore the family of 
rotating minimum energy packets is a distinguished subfamily of all Gaussian states in 
two spatial dimensions, which deserves the detailed analysis. 

The plan of the paper is as follows. Section [2] is devoted to general properties 
of two-dimensional Gaussian packets. The problem of energy minimization under the 
constraint of the fixed mean angular momentum is solved in section [3j The statistical 
properties (such as squeezing, in particular) of the extremal packets and their evolution 
in time are considered in that section as well. In section |4] we calculate variances of 
the energy and angular momentum. The expansion coefficients over the Laguerre- 
Gauss eigenstates (j2| are found in section [5j The difference between ‘co-rotating’ and 
‘anti-rotating’ packets becomes especially clear in quite different expressions for these 
coefficients. In particular, we show that ‘co-rotating’ packets with specific values of 
the ‘internal’ and ‘external’ angular momenta possess more narrow distributions in the 
Fock space than the Poissonian distribution. Section [6] is devoted to the minimal energy 
Gaussian packets for a charged oscillator and a charged free particle in a homogeneous 
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magnetic field. In section [7] we study the peculiarities of the rotating Gaussian packets 
describing a free quantum particle. 


2. Gaussian packets in two dimensions 


We consider normalized Gaussian packets in two dimensions 

ip(x,y) = IV exp [—y (ax 2 + bxy + cy 2 ) + Fx + Gy] ( 6 ) 

where y is a constant scale factor. It is convenient to choose y = Mco/h in the case of 
isotropic harmonic oscillator, but in other cases this factor can take different values. In 


polar coordinates x = r cos( 95 ) and y = rsin(</?) the packet has the form 

- 0 (r, ip) = N exp | — ^r 2 [a + c + (a — c) cos( 2 <^) + &sin( 2 <^)] 

+ r [.Fcos(^) + Gsin(<£>)] j. (7) 

It is useful to separate real and imaginary parts of five complex coefficients 
a = a/2 + i\ai b — f3 + ip, c = y /2 + iy c , F — F\ + iF 2 , G = G\ + iG 2 . ( 8 ) 

Then the probability density has the form 

\^ 2 \(x,y) = \N\ 2 exp[-y(ax 2 + 2/3xy + ^y 2 )], x = x - x 0 , y = y - y 0 , (9) 

where 

\N\' 2 = \N\' 2 exp[y (ax 2 0 + 2/3x 0 y 0 + 'yy/ l )] = yy/A/n, A = ay - f3 2 , (10) 

xq = ('jF 1 -pGJ/iyA), y 0 = (aG 1 - f3F 1 ) /(/iA). (11) 


Obviously a, 7 > 0, whereas all other real parameters can assume any sign, obeying the 
only restriction (3 2 < ay. Parameters x$ and y$ are coordinates of the center of packet. 
They coincide with mean values of coordinates: Xo = (x), yo = (y)- 

The lines of constant relative probability density, defined by the equation 
|^(x,y)| J = exp(—n)|' 0 (xo, yo)| 2 ) are ellipses, whose major/minor semi-axes a±, 
eccentricity e and area Y are given by the formulas 

2v o 2 R TUJ „ n -—-— 


a ± = 


e 2 = 


Y = 


R = y/(a - y) 2 + 4/I 2 . 


a + 7 =F 7? a + y + i? y^A 

The angle 9 between the directions of the major/minor axes and the coordinate axes 
can be found from the equation 


tan(2$) = 2/5 / (y — a). (13) 

Mean values of momenta p x0 = (p x ) and p y0 = (p y ) are 

p x0 = h [F 2 - y (2 Xax 0 + py 0 )] , p y o = h [G 2 - y (2 XcVo + px 0 )\ . (14) 

The probability current density vector has the components 

OV, Jy) = ~v{py + 2x a x],G 2 - n\px + 2\cy]) ■ 


(15) 
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Let us introduce the notation AB = (AB + BA )/2 — (A)(B) for the symmetrical 
covariances of operators A and B. The following expressions hold for covariances of the 
coordinates and momenta operators: 


x 2 =x/(2pA), y 2 = a/(2pA), xy=-fi/(2pA), 

Pi = pH 2 [7 (a 2 + ixl) + a (p 2 - I) 2 ) - 4/3px 0 ] /(2A), 
Pj = pfi 2 [a (l 2 + ixl) + 7 (p 2 - /3 2 ) - 4,9pxJ /(2A), 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


£Pi = R(^P - 27 x„) /(2A), t/p„ =h(fip- 2aXc) /(2A). 


P.P S = pfi 2 p (A - p 2 - 4%aXc) + 2p (ax? + -iXa)} /(2A), 
xp- y = h(2Pxe - Pi) /(2A), gp; = fi (2/3x« - pa) /(2A). 


The mean value of the angular momentum operator L 2 = — yp x can be written 

as (L z ) = hC = h(C c + £i), where the ‘classical’ (related to the motion of the packet 
center) and ‘intrinsic’ (related to the quantum fluctuations) parts are given by the 
following expressions: 


( 22 ) 

(23) 


klCs C X 0 p y0 U()[>x() 1 

Cl = (xpij - Wi) /h = W (xc - Xa) + p{a - 7 )] /(2A). 


It is important that the values of C c and Ci are totally independent for all Gaussian 
packets. Moreover, only packets exhibiting some asymmetry in their shapes (a 7 ^ c) can 
possess a nonzero intrinsic mean angular momentum. Our goal is to find families of the 
‘best’ packets satisfying some additional requirements. 

3. Minimal energy packets of the harmonic isotropic oscillator 

Let us suppose that (J 6 | is the wave function of a particle with mass M moving in the 
isotropic harmonic potential Mix 2 ( x 2 + y 2 ) /2. Then it is convenient to choose the scale 
factor as p = Mix/h. The mean energy £, as well as the mean angular momentum, is 
the sum of two independent terms: £ = £ c + where 



(24) 


and 



Using equations (16)-(18) we can write 


Ci — [(« + 7 ) (l + A + p 2 ) +4 ( 7 X 2 + aXc ) — 4 (3p (x a + Xc)] • (25) 


Obviously the quantities £ c and £ % are totally independent. 
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The minimization of 


he ‘classical’ energy £ c with the fixed value of ‘classical’ 
angular momentum C c ( 22 ) can be done easily: £!f. mri) (C c ) = Jvjj\C c \. The minimizing 


trajectories of the packet center are circles: 

x 2 0 + yl = R 2 = \C c \/n , Pxo = -X c Muy 0 , p y0 = X c Mujx 0 , 
where coefficient X r = ±1 determines the direction of rotation: C r = AJ£J 


(26) 


To find the minimal mean value of ‘intrinsic’ energy £ % for the fixed mean value of 


‘intrinsic’ angular momentum (23) we introduce the new set of parameters 
a + 7 = 2 g, a - 7 = 2f, Xa ~ Xc = 2y, x.a + Xc = 22, 
ot = 9 + f, 7 = 9 ~€, Xa = z + x, Xc = z-x 


and solve equation (23) with respect to y: 
X=(p£-£iA)/(2p). 


Putting (2T)-(29) into (25) we can write £i = £j + E 2l where 

a a' 


Ei(g,v) = ^ kjjjg 


1+ A + 


Eo = -huj 


4 z 2 grj 2 

T + ^a 


p 


V 


££i A 2z/3 


V 


9 


2 z-2 , o2 a 2 2 

V = Z + P , A = g - 9 . 
The minimum of E 2 is obviously achieved for 
2 = 0, p = ^CiA/rj 1 . 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


Therefore the minimum of £i coincides with the minimum of function E\ (30), which is 

(34) 

(36) 


if = 


achieved at (see Appendix A for details) 
9 = 1, 

so that 

^ min) = Mi + IAI)- 



1 Ci 

1 

1 + 1 

Al 


Consequently, the total minimal mean energy equals 

£min — kbJ (1 + \Ci\ + |£ c |) . 


(36) 


If the signs of C, and C c coincide, then £ min = fvjj{ 1 + |£|) in accordance with Q. But 
the mean energy can be much bigger in the case of opposite directions of ‘internal’ and 
‘external’ rotations. 

We see that the minimizing states are degenerate (this is not surprising for the 
isotropic oscillator), since the same values of energy and angular momentum (35) are 
achieved for the Gaussian packets with the following real coefficients (we assume 77 > 0 
in all formulas below): 

a = 1 + ij cos(w), 7 = 1 — 7j cos(rt), y = — Xrjsm(u)/2, 


(37) 
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/3 — r}sm(u), p = Xr]cos(u), 


( 38 ) 


Here A = ±1 is responsible for the sign of mean intrinsic angular momentum (£,; > 0 
for A = 1 and C t < 0 for A = — 1) and u is an arbitrary phase. An additional obvious 
degeneracy is related to the choice of parameters x 0 and y 0 . Complex coefficients of 
function (J6J) have the form (with xq = Rcos(v) and yo = R sin(u)) 


a — - [1 + rj exp(— i\u)] , c = - [1 — r/exp(— iXu )], b = iXp exp(-iXu), (39) 

F = pR {exp (— iX c v) + ?/exp [iX(v — «)]}, (40) 

G = ipR{ A c exp (—iX c v) + A?/exp [iX(v — w)]} . (41) 

Combining equation ([7]) with ( JlTj ) , (14), ( p6| ) and ( p9| we obtain the specihc form 
of the minimal energy Gaussian packets in the polar coordinates: 

ipmin(r,<p) = vVMl - h 2 ) 1/4ex P 


^r 2 [1 + r]exp(2iXip — iXu )] 

$ 

+pRr (exp [iX c (ip - v)] + r/exp[iX((p + v — u)]) - — 


where 


$ = \C c \ [1 + i]cos(u — 2v)\ 


(42) 


(43) 


3.1. The second-order statistical moments and squeezing coefficients 

The covariances of coordinates and momenta in the minimum energy Gaussian packets 
have the following values: 

—77 fli 1 T I T-i I) . r \\ —o h( 1 T \k~'i I)/-, / \\ / a a \ 

x 2 = ——-(1 — 77 cos (u)) , y 2 = —-(1 + r) cos(n)), (44) 


2Mlo 


2Mu> 


pi = 1muH(1 + |£j|) (1 + ijcos(u)), pj = fi(l + |£j|) (1 - tjcos(u)), 
xfa = 7(1 + £,l)Apsin(u), ypy = - 7(1 + £,|)Apsin(u), 
xy = - -ffff L sin(-ti), PxPi = fi(l + |£i|)psin(H), 
xp V = fci [1 - cos(«)/p], yp* = - iliZj [1 + cos(m)/p] . 


(45) 

(46) 

(47) 

(48) 


We see that the partial intrinsic energies £f and £\ coincide, as well as the 
Robertson-Schroclinger uncertainty products U x and U yi where U x = x 2 — (xp x ) 2 : 

u, = u, = L(l + |£,|). (49) 

Let us define the correlation coefficient between variables g and / as r g f = gf / yjg 2 f 2 . 
Then 

rj sin(u) 


r PxPy r xy — ^ V Xp x X r yp y 


yjl — if COS 2 {u) 


(50) 
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The case of u = 0 (pure imaginary coefficient b and maximally different real coefficients 
a and c) corresponds to the absence of correlations in four pairs of variables shown 
above. Another extreme case of u = 7t/ 2 (real coefficient b and complex conjugate 

coefficients a = c*) corresponds to the maximal correlation coefficients in the same 

pairs: \r max \ = 77 = • v /|£|/(l + |£|). 

The best characteristics of squeezing in the x-mode is the invariant squeezing 
,57| (or ‘principal squeezing’ [58] ) 

s x = 2 U x - yjsi - U x j , £ = £/(hcu), U = U/ti 2 . (51) 

For the states under study we have 

S x = S y = < 1, (52) 

so that the states are squeezed, although the concrete ‘directions’ of squeezing can 

be different. For example, for u = 0 we see squeezing in the ^-coordinate and p y - 
momentum. However, the maximal degree of squeezing cannot exceed 50%, since r f < 1. 


coefficient 56 


3.2. Time evolution of packets 


The origin of the degeneracy of minimizing states with respect to arbitrary phases u 
and v becomes clear if one considers the time evolution of these states. It is given by 
the integral 

^(x,y,t)= j G(x,y 1 x',y';t)f(x\y'- I 0)dx'dy', 

where the propagator G(x, y\ x', y’\t) for the two-dimensional isotropic harmonic 
oscillator (calculated for the first time by Kenuard i) reads (here r = (x, y)) 


G(r,r';t) = 




exp 


ip 


[cos(cut) (r 2 + r' 2 ) — 2rr'] > . 


(53) 


2ni sin(cuf) ( 2 sin (cut) 

Performing the integration with initial function ([6]), parametrized as in (39)-(41) with 
initial phases Uo and vq, one can hnd that function ip[x,y\t) has the same form, with 
the only difference that uq and Vo should be replaced by time dependent phases 

u(t) = u 0 + 2A cot, v(t) = v 0 + X c cot. (54) 


This means that the lines of constant probability density are ellipses rotating 
around the central points (x 0 (t), y 0 (t)) with the angular velocity 2A(U without changing 
their shapes (like the ellipses of constant quasiprobability in the phase plane xp of 
one-dimensional harmonic oscillator). The minor axis is inclined by angle u(t )/2 = 
Uq/2 + Xcot with respect to x-axis. The major and minor axes of the ellipse are 
proportional to ( 1 = 1 = 77 ) (with equal scaling factors), and the ellipse eccentricity 
equals £ = [2y/ (1 + g)] 1 ^ 2 . 
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3.3. Verification of universal invariants 

If one combines all coordinate and momentum (co)variances in the symmetric matrix 
Q , then the quantities V m defined according to the expansion 

2 n 

det (Q- 7 £) = ]Tx> m7 m (55) 

m =0 

do not depend on time for any quadratic Hamiltonian. Here £ is the antisymmetric 
matrix constructed from the c-number commutators between coordinates and momenta 
operators, 7 is an arbitrary auxiliary parameter and 2 n the dimensionality of matrices 
Q and £. Quantities V m were named universal quantum invariants in [9l[59 |gT] . They 
are also called as symplectic invariants |62| (or characteristic invariants 163]), especially 
in the quantum information literature [I0[ [l3j[l4|[l6 |. In the one-dimensional case (2x2 
matrices Q and £) the only nontrivial invariant det(Q) coincides with the Robertson- 
Schrodinger uncertainty product U x , and its value for minimal energy packets is given 


by equation (49). In the two-dimensional case (4x4 matrices Q and £) there exist 
two invariants. One of them is again Vq = det(Q). One can verify that the set of 


(co)variances (44)-(48) results in the time-independent value of T> 0 . Moreover, this value 


does not depend on the angular momentum: V 0 = h 4 / 16. Actually, this is the common 
value for all pure Gaussian states, since such states minimize the generalized uncertainty 


relation T) 0 > //'/l(i 19,59, 63,64|. The second invariant of the two-dimensional systems 


V 2 /Tt = ( yp y f + ( xp x f + 2 xp y yp x - 2 xy p x p y 


P 2 X x 2 


P 2 y V 2 


(56) 


also does not depend on the angular momentum for the set (44)-(48): V 2 = — h 4 /2. 
Therefore the combination T>o J r'D 2 /4: + h 4 /lQ equals zero, which is the minimal possible 
value according to another generalized uncertainty relation [9 59 . 

In the current quantum information studies the important quantities are so called 
‘symplectic eigenvalues’ of the covariance matrix, defined as eigenvalues of matrix £ _1 Q. 
It is known that these eigenvalues consist of n pairs («i, — ACi)... (k u , — K n ). Their 
connection with the universal invariants is seen from the identity 

n 

det (Q - 7 E) = det(E)det (E^Q - 7 / 2 n) = det(E) JJ ( 7 2 - k 2 ) , (57) 

3 = 1 

where I 2n is the 2 n x 2 n identity matrix. We have det(E) = h 4 in the two dimensional 
case involved (n = 2). Consequently T >0 = and V 2 = —h 4 + k |). This means 


that the minimum energy states with covariances given by equations (44)-(48) possess 


the minimal possible symplectic eigenvalues \k,\ = \k 2 \ = 1 / 2 , which do not depend on 
the mean angular momentum value C,. 


4. Energy and angular momentum fluctuations 

It is interesting to know the energy and angular momentum variances <Je — {H 2 ) — ( H ) 2 
and ctl = (L 2 ) — ( L ) 2 . One approach is to calculate the fourth order moments of 
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coordinates and momenta. This can be done relatively easy for the Gaussian states, 
because their Wigner functions W (x,y,p x ,p y ) are also Gaussian, so that one can 
use classical formulas (with some modifications due to the non-commutativity of the 
coordinate and momentum operators) for average values of the Gauss distributions 
(see, e.g. [9]). Using the representation of quadrature operators as sums of average and 
fluctuating parts, e.g., x = xq + x, p x = p x o + p x and so on, we can write 

(L 2 ) = (x 0 pyo - yoPxof + Zxoyo [Mu\ c (ypy ~ W^) + (Mufxy - p^Py] 

+ 2Mu)\ c [xpy ( 2x 2 0 + Vo) (%o + ®o)] 

+ {x 2 p 2 y + y 2 p 2 x -xp y yp x -yp x xp y ). (58) 


We take into account that mean values of products of any three operators marked with 
tildes are equal to zero for Gaussian states. The average values appearing in the first 


and second lines of (58) are given by formulas (44)-(48). To calculate the fourth-order 


central moments contained in the last line, we use the known formula for Gaussian 


states, connecting mean values of symmetrical (or Wigner-Weyl) products 65 of four 


operators A , B , C and D (with zero mean values) and sums of pair products of their 
covariances [9] 

(ABCD) W = / W (x,y,p x ,p y ) ABCD dxdydp x dp y /(2iTh) 2 


= AB ■ CD + AC ■ BD + AD ■ BC. 


(59) 


Here A, B, C, D can be any of variables x , y,Px,Py■ The meaning of symbol ( ABCD)w 
is the following: this is the quantum mechanical mean value of the sum of 4! = 24 
products of operators A, B , (7, D taken in all possible orders, divided by the number 
of terms. Mean values of concrete products of operators in predefined orders can be 
expressed in terms of symmetrical mean values with the aid of commutation relations. 
In our case the following relations are useful (they are valid for Gaussian states): 

(x 2 p 2 y ) = 2 {xpj> 2 + x 2 • p 2 , 

(xPyVPx + VPxXPy) = 2(xyp x p y )w + h 2 /2 = 2 (xy ■ p.pp^ + • ypy + Wy • Wc) + ft 2 / 2 - 


After some algebra one can arrive at the following expression for the angular momentum 
variance in terms of the ‘external’ and ‘intrinsic’ mean values C r and £,•: 


f^i/ft 2 — \C c \ + 2 \£i\ (1 + | Gi |) + (1 + AA C ) \C c \ \Ci\ — vWO+M cos(2 w) 
where 


(60) 


w = X(v — u/2) = wo + (AA C — 1) ojt. (61) 

We see that the result depends on the product AA C = ±1, which is positive in the case 
of ‘co-rotation’ of the packet center and ellipse axes and negative for ‘anti-rotating’ 
packets. The phase difference w does not influence the angular momentum variance (as 
well as its mean value) in the ‘anti-rotating’ case: 

or/ft 2 = |U C | + 2 \Ci\ (1 + |£|), AA C = 


- 1 . 


( 62 ) 
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But this phase is important in the case of ‘co-rotation’ (let us assume that £, > 0): 
cr L /h 2 = £ + £,;(! + 2C) — 2£ c y / £7[rT^Ci) cos(2w), AA C =+1, C = C i J rC c . (63) 


The peculiarity of this case is analyzed in section 5.2 


The variance of energy can be calculated in the same manner. Here we need the 
formulas 

x 4 = 3 ( x 2 


x 2 y 2 = 2 ( xy ) + x 2 ■ y 2 , 

(.x 2 pl + p 2 x x 2 ) = 2 x 2 p 2 — h 2 = 4 (xp^) 2 + 2x 2 ■ p 2 — K 2 . 

The final result coincides exactly with ( |60| : 

oe!( hu)~ = \C c \ + 2 \Ci\ (1 + |£j|) + (1 + AA C ) \C c \ \k-'i\ ~ cos(2u;) (64) 


5. Rotating Gaussian packets as superpositions of the Laguerre—Gauss 
energy eigenstates 

It is interesting to find the coefficients of expansion 

V’min ^ ^ C Ur ,mGn r ni 


(65) 


over the Laguerre-Gauss basis ([2]). The simplest formulas correspond to the case 
Ci = y 2 = 0, i.e., the absence of the intrinsic rotation. The state V’mm hr this case 
is nothing but the two-dimensional coherent state with coherent parameters in x and y 


directions a y = ±ia x . Expanding (42) in the Taylor series with respect to the radial 
variable r 




(coh)i 


r,<p) = a/- exp (- (V - 


I Cr. 


£ ^fr exp ^ ikXc ^ ^ ( 66 ) 


k =0 


we see that n r = 0 and the distribution over the energy eigenstates is Poissonian, as 
one may expect for the coherent states. The coefficients are nonzero for the angular 
momentum eigenstates with the same sign of m only: 


Co,k\c ~ 


| r | fc / 2 

Vk\ 


exp (— \C c \ /2 — ik\ c v ), k > 0. 


(67) 


Another simple case corresponds to packets with the fixed center at origin (£ c = 0). 
Then only even azimuthal quantum numbers m = 2kX enter the expansion 

oo i £ 

Vw(A (p) = x/zV^ 1 - r, 2 ) 1/A exp ^ | r \exp [(ik\(2(p - u )]. (68) 

k =0 

Such a structure can be explained by the ‘two-photon’ nature of squeezed vacuum states. 
Comparing (68) with (J2]) we see that nonzero coefficients in expansion (65) are 

— ih\m I Q).2fc+2 

Co,2fcA — 


;-l)‘(l - 


Co,2k 


9 2k + 1 

= v 2 ——-. (69) 


2k+ 2 
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One can notice that coefficients (69) coincide with that of the expansion of the vacuum 


squeezed state of the one-dimensional harmonic oscillator over the Fock basis: 

10 = exp (i [Co' 2 - c-a 2 ]) |0> = (1 - ICD 1/4 E ^r0”|2m> 

k ' rn= 0 


In this case of pure ‘intrinsic’ rotation the vacuum state gives the maximal contribution, 
although the distribution becomes rather flat for states with high mean angular 
momentum (if 1 — rj 2 <C 1 and k 1). 

5.1. Co-rotating packets 

In the general case of £j ^ 0 and C c ^ 0 with equal signs of these mean values (A = A c ) 


we can express function (42) as 


= y^(l - r/ 2 ) 1/4 exp 


h 2 $ 

—r - 

2 2 


£ Miff, (£()/“*•( 70) 


k =0 


k\ 


where Hk(z) is the Hermite polynomial [54] and 
A = e-“”/ 2 VmA B = [rje*" + e - *"] V|£ c |/(2ij), w = \(v - m/2). 


(71) 


Consequently nonzero coefficients of expansion (|65|) have the form 
Co.fcA = 


(1 -’ ,2)1/4 (e-“V2)‘ /2 






7e™ + 


]v/H7M 


(72) 


x exp {— \C c \ [1 + g cos(2w)] /2} , k > 0. 

Note that phase w does not depend on time for coinciding directions of the ‘internal’ 
and ‘external’ rotations. 

To understand better the influence of phase w on the behavior of the probabilities 
Pk = |Co,fcA | 2 , let us consider first the special case of w = 0 and C > 0: 


Pk = 


(1 — rj 2 ) 1 / 2 ^ 

2 k k\ 


hi (h+ijyz; 


exp [-£ c (l + i?)] • 


(73) 


For big ‘internal’ angular momenta C L > 1 we consider t] = 1 — e with e < 1, so that 
Ci ~ (2e)^ 1 . Then equation (73) can be simplified if Ci C c . Using the approximate 


formula (1 — e) k ~ exp(— ke) we get the expression 

Pk ~ exp t _2£c _ k / 

where contributions of ‘internal’ and ‘external’ rotations are factorized. In the most 
interesting region k ~ C 1, the argument of the Hermite polynomial is much smaller 


(74) 


than its index. Then the known asymptotics of the Hermite polynomials 54 55 together 


with the Stirling formula for factorials lead to the most simple expression if £; > £ c > 1: 

2 exp {-k/C) _ 2 ( in r TFT , ^ /o\ ^ 


Pk 


vS cos 2 (v / 2£ c (2fc + l) -W2). 


We see that probabilities pk rapidly oscillate around some slowly decaying average 
distribution. Replacing cos 2 (...) by its average value 1/2 and integrating over k from 0 
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to oo (i.e., using the simplest form of the Euler-Poisson summation formula) we arrive 
at the correct normalization = 1- This shows the reliability of the approximate 

formula (75). 

Now let us suppose that w = 7t/ 2, £,. 1 and £ c -C £*. Then the argument of 

the Hermite polynomial in equation (72) equals approximately —isJCJ (8£ 2 ), i.e, it is 
close to zero. The argument of the exponential function in (72) equals approximately 
—£ c / (2£j), so that it is close to zero as well. Consequently, the distribution p k is close 
to ( J69| ) in this case. 

ft can be interesting to calculate the mean value of the angular momentum using 
coefficients (72). This can be done with the aid of generating function (dependent on 
an auxiliary variable z) 

OO 

G ( z ) - i c °> 


i,fcA 


v = 


k =0 


1 — I) 2 

1 — z 2 rj 2 


exp 


\Cc\ (z - 1 ) 


1 — Z 2 T] 2 


[l — zrj 2 + 77 ( 1 — z) cos(2w)] 


(76) 


We used here the known Mehler formula 

(C/2) 1 


£ 

k =0 


k\ 


-H k (x)H k (y) = (1 - C 2 ) 


- 1/2 


exp 


2xy( - (x 2 + y 2 )( 2 

i-C 2 


The value G(l) = 1 confirms the correct normalization of coefficients ( |72| ). The mean 
value of the angular momentum can be calculated as k |c 0 ,/eA| = hG'{z)\ z= i- The 

result (for £,; > 0 and £ c > 0) is h [£ c + if /(1 — if)} = K (£ c + £j) = TiC independently 
of the phase w. The variance of the angular momentum a k = (L 2 ) — (L z ) 2 is determined 
by the first and second derivatives of G(z) at z — 1: Ui/h 2 = G"( 1) + G'(l) — [G'(l)] 2 . 
The result coincides with equation (63). 


5.2. Competition of ‘external’ and ‘ intrinsic’ rotations in the co-rotating case 

If L % — 0 (the rotation of coherent state circular-shape packet along a circle) then 
we have the typical result for the poissonian distribution ctl = h 2 C. In the case of 
C c = 0 we have the typical formula for fluctuations in the squeezed vacuum state 
(Tl = 2h 2 |£|(l + |£|). But it is remarkable (and perhaps unexpected) that fluctuations 
of the angular momentum can be smaller than that in the Poissonian distribution if both 
quantities £, and £ c are different from zero (some analog of sub-Poissonian statistics). 
Indeed, for £j —* 0 we have ai/h 2 ~ £ [l — 2v / £icos(2w)] < £ if cos(2ic) > 0. The 
maximal ‘squeezing’ of the distribution |Co,fcA| is observed for w = 0, when the minor 
axis of the coordinate probability density ellipse is directed along the radius connecting 
the center of packet with the origin (i.e., the ellipse is ‘squeezed’ in the radial direction): 
see figure [lj 

The minimum of ap with the fixed value of mean angular momentum hC and w = 0 
is achieved for £* satisfying the equation 

C = y'cj] - £,] (1 + SC, + 8 Cf) + SC, + 12 C\ + 8 C\. 


(77) 
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Figure 1. The orientation of the ellipse of constant probability density with respect 
to the circle which the center of packet moves along in the case of relative angle w = 0. 
This relative configuration is maintained in time in the co-rotating case CiC c > 0. 


Combining (63) and (77) we obtain the minimal variance of the total angular momentum 


as function of the ‘intrinsic’ momentum £p 


cr. 


nn /hr — 4 Ci (1 + Ci ) + (1 + 2 Ci) \JCi (1 + Cj). 


(78) 


We see that for £ t <C 1 the minimum is achieved if £, ~ £ 2 <C 1. Then (7™ m /h 2 ~ £—£?. 
The optimal ellipse eccentricity in this case equals £* ~ \/2£ <C 1. Actually, this 
happens for very small values of £ % . Increasing £, and respectively £, one can make the 
effect of diminishing the variance (compared with the case of rotated coherent packet) 
more significant. For example, taking £ t = 1/8 we obtain £ = 13/8, a™ in /h 2 = 33/32 
and e* = l/\/2. For £ % = 1/3 we have £ = 19/3, a™ in /h 2 = 26/9 and £* = ^2/3. 


An approximate solution to equation (77) for £ % » 1 is £ « 16£f. Then 
L I n ~ 6£ 2 ~ (3/2)(£/2) 2 / 3 . Consequently the relative minimal width of the angular 
momentum distribution 


aT n /h 2 » 


yja™ in /(h £) « ^/3/2V3(2£)- 2 / 3 (79) 

can be significantly smaller than the Poissonian width £~ 1//2 in this asymptotical case, in 
spite of the small relative weight of the ‘intrinsic’ angular momentum £i/£ ~ (4£) -2 / 3 . 

5.3. 1 Anti-rotating ’ packets 


To find the coefficients of expansion (|65|) in the case of A = — A c we represent the 

(80) 

*(r, v) = £ e “ v £' ”'' 

where 


exponential function in the right-hand side of (J42J) as the product of two sums: 
■0 min(r,<p ) = a//Vtt( 1 - ?? 2 ) 1/4 exp (-/xr 2 /2 - <h/2) T(r,<^), 

~ M)* 


fc =0 




3=0 


3 ! 


(81) 


Si = (|£ob/2) 1/2 e‘”. 

Combining the terms with the same powers of e %Xip we can write 

oo oo / 2 n \i 

*(r, U = £ e ln»A(r— u/ 2 ) (,, r U/2)’" /2 V )// 1 (S,) 


(82) 


ra=0 


pj j' ! 0 +i) ! 


£■ 

m=l 




/ 2 , ,xm/2 Ip-'-Dp „ r R 

£r?7VTT7T B ‘(B. 


(/ir 2 5i Nfc 


jP /d(m + fc)! _ 


(83) 
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Using the formula [55] (we consider integral positive values of m here) 


x k = ST' (~ l) n k\(k + m)\ (m) 


^ (n + m)!(fr-n)! 


(84) 


we continue as 


V(r,<p) = Y e imKip - u/2) 


-B.r r , mW 


ra =0 

oo 


_n V ^ r j= o 


m+n+ji (-Bl 


n=0 

oo 


e —im\(tp—v) ( /ir 2 |£ c |) m/2 




(85) 


m=l 


-^(n + m)! n v '/ji 
n=0 ' ' fc=0 


Comparing ( |80| and (85) with (|2]) and using formula 5.12.1.3 from |66 

°° y.fc 

E ry-Hn+fc^) = exp (2xf - f 2 ) iP n (x - t) 


fc=0 


( 86 ) 


we arrive at the following expressions for coefficients c nr)m \ 

_ (l-^) 1 / 4 (-B 1 )’ 1 - ^_ K .l /2n 

°n r ,Am / pr-j p-r ° J - y n r ,m’> 

xjn r \ ( n r + |m|)! 
where <p = sin(2u>) \C c \ y/2 and 

_/ (r)e- iXu /2) m/2 H m+n M, rn > 0 

(\C c \e 2iXv ) lml/2 H nr (0), m< 0 

We see that only even radial quantum numbers n r give nonzero contributions for m < 0, 
whereas the parity of n r must coincide with the parity of m for m > 0. Remember that 


Dn r: rn 


(87) 


( 88 ) 


H 2k { 0) = (-1) 


k m 


k\ 


> Vii;+l(0) — 0. 


(89) 


2 

Probabilities \c nr ,\m\ do not depend on phases v and u. They satisfy the normalization 
condition, which can be written as (we use here the notation x = r]/2 and y — \C c \ to 
simplify the formula) 


OO OO 


EE 


—' ^ n\(m + n)\ 

=o m= 1 v ’ 


V) (x m [H m+n (0)} 2 + y m [H n (0)} 2 ) = 


\/l — 4x 


(90) 


This identity is not obvious at first glance, but it can be proven after some algebra, ft is 
interesting that in spite of strong ‘entanglement’ between the ‘intrinsic’ and ‘external’ 
rotations (characterized by parameters y and C c ) in the formulas for coefficients c nrt \ m , 
these rotations are totally disentangled in formulas for the mean angular momentum, 
mean energy and their variances, as was shown in the preceding sections. 
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6. Isotropic charged oscillator and charged particle in a homogeneous 
magnetic field 


The stationary Schrodinger equation for the two-dimensional isotropic oscillator in a 
homogeneous magnetic field B , described by the Hamiltonian (in the circular gauge of 
the vector potential) 


H = 


1 -s 2 | M -2^1 f 

—ttP T —oj r — u )£, 1 j z , 
2 AT 2 ’ 


UJ L = 


eB 
2Me’ 


~ 2 2,2 
lu — ca T c 


was solved in polar coordinates by Fock 67 
(r, <p) = 


fa 


fin r 




A 4 „2 


(91) 


(92) 


if nrm = hu (1 + |m| + 2n r ) — hupm- 


(93) 


The solution (92) differs from (J2| by the change fi —* jl = MCu/h. The extremal Gaussian 
states with the fixed values of L % and C c are given by formula (42) with parameter /j 
replaced by jl. Their mean energy equals 


£ = hio (1 + \Cj\ + |£ c |) — hu>L (Ci + jC c ) 


(94) 


The minimal energy is achieved for co-rotating packets with equal signs of C and ujl- 


£min{£) = kuu + h(u - \u L \) |£|, > 0. 


(95) 


Such packets are superpositions of energy eigenstates with n r = 0 and coefficients (72). 

In the case of charged particle in a homogeneous magnetic field {uj = 0, Cj = \ojl\) 
we get (assuming uj l > 0) 

£ = frwL [1 + \Ci\ (1 — A) + \C c \ (1 — A c )]. (96) 


The absolute minimum £ min = H\ujl\ is achieved for all co-rotating packets with 
Cujl > 0. Of course this is explained by the well known infinite degeneracy of energy 
eigenstates in this special case. 

The energy variance can be calculated in the same way as in section [IJ The result 
is (for ul > 0) 


a E /(hw L ) 2 


2 (1 - A c ) (1 - A) \C c \ [|A| - VlA|(l + |^|) cos(2 W ) 
+ 2 \C c \ (1 — A c ) + 4 \Cj\ (1 + |/lj|) (1 — A). 


(97) 


The variance equals zero for all packets whose directions of ‘internal’ and ‘external’ 
rotations coincide with the direction of the Larrnor rotation: A = A c = 1. The relative 
phase w is important if only A = A c = —1 (packets performing ‘co-rotation’ in the 
direction opposite to the Larmor rotation). For ‘anti-rotating’ packets, the energy 
variance equals either 4 \C c \ (if A c = —1) or 8 \Ci\ (1 + |£j|) (if A = —1). 

Applying the propagator |2| (with jij J = Mujl/K) 


G(r,r']t) 


Rl 

2m smjuiLt) 


exp 


ifJ-L 

2 


cot(c oi,t) (r — r 7 ) 2 


2 to/ 



(98) 
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to the initial function ([6]), parametrized as in (39)-(41) with initial phases Uq and Vo, 
one can see again that function ip(x,y\t) maintains its form, provided w 0 and Vq are 
replaced by time dependent phases 

u(t) = uq + 2u L t(\ - 1), v(t) — v 0 + u L t (A c — 1). (99) 

Thus we see again that all packets with A = A c = 1 do not rotate at all, although they 
can possess arbitrary values of ‘external’ and ‘internal’ angular momenta (but the same 
minimal possible total energy hi ol)- 


7. Free particle 

In the case of free particle, the center of packet moves along a straight line, so that 
there is no ‘external’ rotation. Therefore it is sufficient to consider the special case of 
homogeneous Gaussian packets with the center fixed at point (0, 0). There is no positive 
minimal energy for a free particle, therefore the results of this section cannot be obtained 
as a limit oo —> 0 of the preceding sections. The qualitative difference is that the shape 
and size of free packets are inevitably deformed during their evolution in time, whereas 
the minimum energy packets studied in the preceding sections maintain their shape and 
size. For this reason we put hereafter /j = 1, using dimensional coefficients a, b and c. 

The well known free particle propagator 2| 

irn , (N 


G(r,r']t) = [2iriht/m] 11/2 


exp 


2 ht 


( 100 ) 


transforms initial packet (J6J) to another Gaussian packet with time dependent coefficients 
a' = (a + irD)/G(T), d = (c + irD)/G(r), b' = b/G(r), (101) 

where 

r = 2ht/m , D = ac — b 2 / 4, G(t) = 1 + ir(a + c) — t 2 D. (102) 

To diminish the number of parameters, let us consider the initial ‘most symmetrical’ 
packet with nonzero mean angular momentum: a = 7 = a 0 , Ac = ~Xa = Ao > 0, /3o > 0 
and po = 0. Then the conserved mean angular momentum is positive: 

2A)Xo 


C = 


(103) 


<4 - K 

The ellipses of constant probability densities have the semi-axes a^_(0) = vj (ap =F A)), 
major semi-axis being inclined at the angle —7 t/ 4 with respect to ir-axis. 

However, the symmetry of coefficients is destroyed in the process of evolution, 
because real and imaginary parts of complex coefficients ( 101[ ) have the following form: 

(104) 

(105) 


a(r = 


01 0 

W: 


(l + t 2 D + - 2rxo) , 70) = (l + t 2 D + + 2txo) , 


F{r) 


do 


P(t) =-—-(l-T 2 D + ) , p(r) = — 


/3 0 a 0 T 


F(r) 

Xa,c( T ) = [TAO (1 - T 2 D + ) 


D_t 


Hr) 

D\r 3 ] /F(t 


A (r) = 


A(0) 

F(t) 


(106) 
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where 


F(t) = 1 + 2t 2 £>_ + D\t‘, D ± = - (ajj ± 4xl =F P 2) 


(107) 


Note that functions G(r) and F(t) introduced above have no relation to the coefficients 
G and F used in section [2] 

Contrary to the harmonic oscillator case, the free particle ellipse of constant relative 
probability not only rotates, but it changes its shape. This is clearly seen from the 
formula describing the ellipse area Y(t) or the probability density at the center of packet: 

Y (j)/Y(Q) = |-0(O, 0; O)/-0(O, 0; r )| 2 = \f~F{j). (108) 

While coefficient D + in equation ( |107 ) is always positive, coefficient _D_ is negative if 
4 %q > ckq + ^o- In such a case the packet shrinks initially, since 1-0(0, 0; 0)| J < |0(O, 0; r )| 2 
for small values of r. Such a behavior is not surprising for correlated (y 0 7 ^ 0) Gaussian 
packets in one dimension, since parameter y, depending on its sign, is responsible for 
the effects of focusing or de-focusing (therefore correlated free Gaussian packets in 
one dimension were named ‘contractive states’ in [68l 691). This is clearly seen from 


formula (15), which shows that the initial probability current is directed to/from the 
center if y ^ 0. A nontrivial effect of quantum shrinking of packets with initial zero 
probability current density was discovered in 70 and generalized in 71 75 . Note that 


packets studied in this connection in two dimensions had the ring-shaped forms with 
zero probability density at the origin 70,72,73 . We see that the 2D Gaussian packets 


with the maximum of probability density at the origin can also shrink, but this time this 
happens because of some competition between focusing and de-focusing in orthogonal 
directions (y a = — y c ), and this happens for sufficiently strong focusing only. 

The minimal value 

4ao (4%o — fio) 


F = 

J- mm. 


H + 4 xl - /3q) 2 


(109) 


is achieved for rA n = —D~/D 2 + . This value can be made as small as desired for big 
enough values of parameter xo (so that the packet can be concentrated in a very small 
region at r = r min ). Note that this can happen even for /3 0 = 0, i.e., for zero mean 
value of the angular momentum. This shows the importance of parameter xo, which 
determines the initial correlation coefficient between the coordinates and conjugated 
momenta. 

If /3o 7 ^ 0 (i.e., C 7 ^ 0), then some rotation of the packet is observed: 

2/5 0 (1 - t 2 D + ) 


tan[ 2 d(r)] = 


( 110 ) 


4txo«o 

The major axis becomes parallel to one of coordinate axes when t 2 D + = 1 (one can 
easily check that this happens after r mm if r min exists). After that instant, angle 6 
changes its sign, going asymptotically to the same absolute value 7 t/ 4 as at t = 0. This 
means that the asymptotic direction of major axis is perpendicular to the initial one. 
The ellipse also changes its shape in the process of evolution: 

Aq 


*± 


,T = 


a 0 (1 + t 2 D + ) ± \J$ {1 - r 2 D + f + 4a§x§r 2 


( 111 ) 
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For 1 we have 


*± 


(r)«4(0) 1 + 


T 


4|/5o 


[|/5o| («o =F |/?o|) 2 ± 4xo (2«o =F |/?o|)] 


Consequently, the major semi-axis always increases, whereas the minor one decreases 
for big enough values of | x' 0 1 • I* 1 particular, for sure it decreases for 4xq > «d + /3q. 

To understand the evolution of the ellipse eccentricity we transform the 
corresponding expression in (12) as follows: 


2 _ i « + 7 
E 2 R 


= 1 + 


4A 


(a + y) 2 _ 


- 1/2 


In turn, using (104) and (105) we have 


4A 


(a + y ) 2 


a; 


A(0)F(r) 
j(l + r*£> + ) s 


A(0) 


cy,r 


1-2 (D, 


D) D 


1 

-T + - 

r 


-2 


Consequently, the eccentricity attains the maximal value 


^ m n 'r' 


21 To 


4|Xo 


I To I + VD+ 21 xo | + \J 4\o + « 2 - P 2 


( 112 ) 


exactly at the time instant To when the ellipse axes become parallel to the coordinate 
ones: TqD + = 1. At this instant we have 


i ± 


fo) = 


2a 0 a±(0) 
«o ± |/3 0 | 


(■*»)• 


FO o) = 


4ccq 


+ 4Tq - Po 


(113) 


VDPJ ' " «6 

Obviously F{tq) > F min . Nonetheless F(t 0 ) still can be smaller than unity for big 
values of |y 0 |- In such a case, the area of ellipses at r = r 0 is smaller than the initial 
one. Function F(t) returns to the initial unity value when r 2 = —2 D_/D\ = 2r 2 wn . 

When r oo, then a±(r) ~ a±(0 )y/D+ r and the eccentricity returns to the 
initial value e 2 (0) = £ 2 (oo) = 21/3 0 1/ («o + |/?oI)- Tde probability density at the center 
decreases asymptotically as |-0(O, 0; r)/-0(O, 0; 0)| 2 ~ ( D + t 2 ) . The change in time of 
the shape and orientation of the ellipses of constant relative probability density is shown 
schematically in figure [2] 


8. Conclusion 

We have studied the properties of two-dimensional Gaussian packets with fixed mean 
values of angular momentum. They depend on mutual directions of the independent 
‘internal’ and ‘external rotations’ in the case of two-dimensional harmonic isotropic 
oscillator. Moreover, the direction of the Larmor rotation is also important in the 
presence of an additional homogeneous magnetic field. The states minimizing the total 
mean energy possess nonzero correlation coefficients between coordinates and conjugated 
momenta. They also show a moderate squeezing of the quadrature components. The 
distribution function over energy eigenstates can exhibit some kind of ‘sub-Poissonian’ 
statistics for ‘co-rotating’ packets satisfying certain conditions between the mean values 
of ‘internal’ and ‘external’ angular momenta. In the case of free particle, packets with 
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Figure 2. The orientation and shape of the ellipse of constant relative probability 
density of the free Gaussian packet with big initial coordinate-momentum correlation 
coefficients at four specific time instants. 1 - the initial instant r = 0. 2 - the instant 
Tmim when the probability density at the center attains the maximal value and the 
ellipse area is minimal. 3 - the instant r 0 , when the eccentricity attains the maximal 
value and the ellipse is rotated by 45° with respect to the initial position. 4 - the 
asymptotical behavior: the eccentricity returns to the initial value, the ellipse is rotated 
by almost 90° and its area grows unlimitedly. 


big enough initial coordinate-momentum correlation coefficients shrink initially. Only 
after some time they start to expand, rotating the directions of their major/minor axes 
of constant probability density ellipses by 90 degrees. Since the shrinking effect can be 
very strong for the packets with big initial coordinate-momentum correlation coefficients, 
such packets, perhaps, could find applications in the sensitive electron microscopy. The 
time T m i n of the ‘maximal shrinking’ can be adjusted to the necessary distance between 
the beam source and target by means of choosing the appropriate longitudinal velocity 
of the real three-dimensional beam along the propagation axis. 
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Appendix A. Proof of solution (34) 


Using the relation A = g 2 — g 2 one can represent the equation dE\/dg = 0 as 
g/A 2 = g 2 C 2 /g 3 , so that g 2 = \Ci\gA (remember that g > 0 and A > 0). Then it 
is easy to obtain the relations 


g 2 = 




A = 


9 


1 + 9 1 £■; I 1 + 9 1AI 

The equation dEi/dg = 0 has the form 

1 1 C 2 A o 2 (C 2 1 , n 

1 + A + ~^ + 2g (?“ Ad -°- 


(A.l) 


(A.2) 


Putting expressions (A.l) into the left-hand side of (A.2) one arrives at the simple 
equation 1 — 1/g 2 = 0, resulting in the solution (34). 
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